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Abstract 

In a sample of 58 million J/ip events collected with the BES II detector, the process 
J/ip — > 77/ c is observed in five decay channels: rj c — ► K + K~ir + ir~ , ir + ir~ ir + ir~ , 
K^KgH^ (with Kg — > 7r + 7r~), (pep (with cp — ► K + K~) and pp. From these signals, 
we determine 

Br(J/^ -► 777c) x Br(r] c -► K^K'^-k') = (1.5 ± 0.2 ± 0.2) x 10" 4 , 
Br(J/ip -> 7%) x Br(r] c -> vr+vr^vr+Tr") = (1.3 ± 0.2 ± 0.4) x 10~ 4 , 
Br(J/tp -► 777c) x Br{r] c -► K^^ttF) = (2.2 ± 0.3 ± 0.5) x 10" 4 , 
Br(J/i) -► 7 77c) x 5r(r/ c -► 00) = (3.3 ± 0.6 ± 0.6) x 10~ 5 and 
Br(J/tp 7 77c) x .Brfr/c -»• pp) = (1.9 ± 0.3 ± 0.3) x 10~ 5 . 

PACS: 13.25.Gv, 14.40.Gx, 13.40.Hq 
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Hadronic decays of the r\ c have been studied by Mark III [1,2], DM2 [3], and 
other experiments [4-7]. However, the branching fractions of the i] c still have 
very large errors in the Particle Data Group (PDG) compilation [8]. More 
recently the branching fractions for B — > rj c K decays and B — * r] c K* have 
been measured by the Belle [9,10] experiment, and their measured branching 
fraction for rj c — > 00 is smaller than the PDG value [8]. 

In a previous paper [11], based on 58 million J/ip events collected in the 
Beijing Spectrometer (BES II) detector at the Beijing Electron- Positron Col- 
lider, we measured the r\ c mass and width using the processes J/-0 — * 7%, 
r] c -> K+K-iT+iT-, 7r + 7r-7r + 7r-, K± K%ir* (with K° s -> 7T+7T-), 00 (with 
— > K + K~) and pp, and obtained m^ c = 2977.5 ± 1.0 (sta) ± 1.2 (sys) 
MeV and T Vc = 17.0 ± 3.7 (sta) ± 7.4 (sys) MeV. In this paper, we report 
measurements of the branching ratios for the same processes. 

BES is a conventional solenoidal magnet detector that is described in de- 
tail in Ref. [12]; BESII is the upgraded version of the BES detector [13]. A 
12-layer vertex chamber (VTC) surrounding the beam pipe provides trigger 
information. A forty-layer main drift chamber (MDC), located radially out- 
side the VTC, provides trajectory and energy loss (dE/dx) information for 
charged tracks over 85% of the total solid angle with a momentum resolution 
of <Jp/p = 0.0178\/1 + p 2 (p in GeV/c) and a dE/dx resolution for hadron 
tracks of ~ 8%. An array of 48 scintillation counters surrounding the MDC 
measures the time-of-flight (TOF) of charged tracks with a resolution of ~ 200 
ps for hadrons. Radially outside the TOF system is a 12 radiation length, lead- 
gas barrel shower counter (BSC). This measures the energies of electrons and 
photons over ~ 80% of the total solid angle with an energy resolution of 
(Te/E = 21%/y/E (E in GeV). Outside the solenoidal coil, which provides a 
0.4 Tesla magnetic field over the tracking volume, is an iron flux return that 
is instrumented with three double layers of counters that identify muons of 
momentum greater than 0.5 GeV/c. 

A Geant3 based Monte Carlo, SIMBES, which simulates the detector response, 
including interactions of secondary particles in the detector material, is used in 
this analysis. Reasonable agreement between data and Monte Carlo simulation 
is observed in various channels tested, including e + e~ — > (-y)e + e~, e + e~ — > 
(7)/x/x, J/tp — > pp, J/4> — > pvr and ip(2S) — > n + n^J/ip, J/ip — > l + l~. 

The event selection criteria for each channel are described in detail in our pre- 
vious paper [11]. Here we repeat only the essential information and emphasize 
those considerations that are unique to the r\ c branching ratio measurement. 

Candidate events are required to have the correct number of charged tracks 
for a given hypothesis. Events are kinematically fitted with four constraints 
(4C) to the hypotheses: J/ip — > 'jK + K~it + it~ , J ftp — > 77r + 7r~7r + 7r~ , J/ip — > 
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7-ft' ± 7r T 7r + 7r~, and J/ip 7pp. A one-constraint (1C) fit is performed for the 
J/ijj — > r f m i SS K + K~ K + K~ hypothesis, where 7 m «ss indicates that this photon 
is not detected. Events with a x 2 l ess than 40.0 for a particular channel are 
selected. 

In order to remove backgrounds from non-radiative decay channels, all se- 
lected events are subjected to (4C) kinematic fits to the hypotheses: J/x/j — > 
K + K^n + n^ , J/tp — > 7r + 7r~7r + 7r - , and J/^ — > _ft' ± 7r T 7r + 7r _ and are required 
to satisfy xV/^ -> K+R-n+n-) > 20.0 (for K+R-n+n-); x 2 (J/*P -" 
7r + 7r _ 7r + 7r _ ) > 10.0 (for 7r + 7r~7r+7r~) and x {J/^ — > K ± n' f 'K + 7r~) > 10.0 
(for fr ± i^57r =F ). For the J/^ — > 7pp channel, we require that the opening an- 
gle of the two charged tracks is smaller than 179°. A detailed Monte Carlo 
simulation shows that these cuts, referred to below as the J/ip veto, do not 
distort the invariant mass distributions around the r\ c signal peak. 

After event selection, the invariant mass spectra for the individual decay modes 
are obtained, as shown in Fig. 1. An unbinned maximum likelihood fit using 
MINUIT [14] is performed to all five channels simultaneously. The fitting 
method is described in detail in our previous paper [11]. 

The branching ratio can be calculated using 

Br = N I g 1 = JL t 

where e is the detection efficiency; N = Nf it /e is the efficiency-corrected num- 
ber of 7] c events obtained directly from the fit and corrected using Br(K® — > 
7T+7T-) and Br((j) -> K + K~) [8] where necessary; and Nj/^ = (57.7 ± 2.72) x 
10 6 [15] is the total number of J/ip events. The numbers of r] c events deter- 
mined from the fit and the corresponding branching ratios, by decay channel, 
are listed in Table 1. 
Table 1 



Number of rj c events and corresponding branching ratios for the individual channels 
(corrected using Br(K® — > tt + tt~) and Br{4> — ► K + K~) [8] where necessary). 



Process 
J/ip -► 7r? c , 


No. of events 
(detected) 


No. of events 
(efficiency-corrected) 


Product of 
branching ratios 


r] c -► K+K--K+-K- 


413 ± 54 


8453 ± 1110 


(1.5 ±0.2 ±0.2) x 10~ 4 


r\ c — > 7r + 7r~7r + 7r~ 


542 ± 75 


7643 ± 1062 


(1.3 ±0.2 ±0.4) x 10" 4 




609 ± 71 


12516 ± 1460 


(2.2 ±0.3 ±0.5) x 10~ 4 


»7c -»■ 4>4> 


357 ± 64 


1922 ± 357 


(3.3 ±0.6 ±0.6) x 10~ 5 


Vc ~^PP 


213 ± 33 


1105 ± 171 


(1.9 ±0.3 ±0.3) x 10" 5 



The main systematic error contributions in measuring the r] c branching ratios 
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Fig. 1. Invariant mass distributions in the r\ c 

"ivr+TT-Tr+Tr-, (c) m K±K o 7T ^, (d) and (e) m, 

the data; the curves are the fit result. 



region (a) m K+K - 7T+n -, (b) 
pp. The histograms correspond to 



originate from uncertainties in the background shape parameterization used, 
differences between different Monte Carlo simulations of the drift chamber 
wire resolution, detection efficiency differences due to uncertainties in i] c de- 
cay sequences into the final state (for i] c — ► 7r + 7r~7r + 7r~, r] c — ► K + K~7t + tt~ 
and 77 c — > K ± Kg7T^), differences in the photon efficiency determined using 
data and that determined from the Monte Carlo simulation, particle identifi- 
cation uncertainties, and the uncertainty in the total number of J /if) events. 
In Fig. 1, second-order polynomials are used to describe the backgrounds. The 
systematic errors due the background shape are studied by using instead lin- 
ear polynomial functions to fit the backgrounds in Fig. 1(b), (d), and (e) and 



5 



third order polynomials to fit the backgrounds in Fig. 1(a) and (c), changing 
the upper fitting bound from 3.05 to 3.07 GeV/c 2 , and removing the J/vp veto 
from the event selection. The relative systematic errors from these sources are 
listed in Table 2. Since the errors are correlated, we choose the largest one as 
the systematic error due to the background shape. 

Table 2 



Relative systematic error caused by background shape. 



Sources 


K+K-TT+1T- 




K±K S ^ 




PP 


background polynomial 


4.4% 


7.6% 




2.5% 


8.3% 


3.2% 


fitting range 


9.4% 


8.4% 




17.2% 


15.5% 


10.6% 


J/ip veto 


1.7% 


26.6% 




10.1% 


17.3% 


15.2% 



The relative systematic errors for the individual channels are summarized 
in Table 3, where the individual contributions are added in quadrature to 
obtain the total relative systematic error. The systematic errors on the product 
branching ratios are given in Table 1. 

Table 3 



Relative systematic error summary. 



Sources 


K+K-1T+7T- 


7T + 7T 7T + 7T 


K±K° s tt^ 


00 


pp 


BG shape 


9.4% 


26.6% 


17.2% 


17.3% 


15.2% 


wire resolution 


10.4% 


17.1% 


13.1% 


2.9% 


4.7% 


7] c decay sequences 


4.5% 


4.5% 


1.0% 






7 efficiency 


2.0% 


2.0% 


2.0% 


2.0% 


2.0% 


particle ID 


2.5% 


2.7% 


2.2% 


2.5% 


1.1% 




4.7% 


4.7% 


4.7% 


4.7% 


4.7% 


Total 


15.8% 


32.5% 


22.3% 


18.4% 


16.7% 



Using the branching fraction Br(J/ip —± 7%) = (1.3 ± 0.4)% [8], the f] c 
branching fractions can be obtained. Table 4 shows the BES results together 
with the PDG [8] and Belle [9,10] values. The BES Br(r] c -> 00) is smaller 
than the current PDG value of (7.1 ± 2.8) x 10~ 3 and is consistent with the 
Belle [10] and DM2 [3] measurements within errors. The branching fractions 
for i] c — > K ± Kg7i T and i] c — > pp are consistent with both the Belle[9] and PDG 
values [8]. The branching fractions for r\ c — > 7r + 7r~7r + 7r~ and r\ c — > K + K^n + n^ 
are consistent with the PDG values [8] within errors. 
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Table 4 

Branching fractions of the rj c (the Belle results of Br(rj c — > ET ± i ; rg7r =F ) and Br(rj c — ► 
pp) are calculated from reference [10]). 



Process 


BES(%) 


PDG02(%) [8] 


Belle(%) 


Br(r] c -> ET+K-vr+vr-) 


1.2 ±0.4 






Br(r] c — > 7T + 7r _ 7r + 7r _ ) 


1.0 ±0.5 


1.2 ±0.4 




Brfe -> i^if^) 


1.7 ±0.7 


±(5.5 ±1.7) 


~ 1.8 




0.25 ±0.09 


0.71 ±0.28 


0.18±g;g| ±0.07 


Br(r] c — > pp) 


0.15 ±0.06 


0.12 ±0.04 


~ 0.14 
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